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I. INTRODUCTION 
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A non-perturbative renormalization of a many-body problem, where non-relativistic bosons living 
on a two dimensional Riemannian manifold interact with each other via the two-body Dirac delta 
- . . potential, is given by the help of the heat kernel defined on the manifold. After this renormalization 

procedure, the resolvent becomes a well-defined operator expressed in terms of an operator (called 
, principal operator) which includes all the information about the spectrum. Then, the ground state 

' energy is found in the mean field approximation and we prove that it grows exponentially with 

(N ■ 

the number of bosons. The renormalization group equation (or Callan-Symanzik equation) for the 
^ , principal operator of the model is derived and the (3 function is exactly calculated for the general 

O ' case, which includes all particle numbers. 

Q 

o 

m 

' Ultraviolet divergences appear not only in quantum field theories [l[ but also in many-body theories and non- 
2 I relativistic quantum mechanical problems in vifhich the interaction has a peculiar singular behavior at short distances 
■ In all these cases, infinities are encountered when we calculate some obscrvablcs (experimentally measured 
quantities) e.g., differential cross section of a scattering process, bound state energy, etc. In order to circumvent 
these divergences, a series of algorithmic steps must be applied, and this whole procedure is called renormalization. 
^ ■ The basic idea of renormalization is first to regularize the infinite integrals by modifying the short distance (or large 

■ momenta) behavior of the interactions for ultraviolet divergences. This can be accomplished in several ways with 
f — . ' the assumption that the theory is valid up to a scale determined by an unknown parameter, called cutoff e (or A in 
T-H J momentum space). According to the modern point of view of renormaUzation Q, a renormalizable theory could be 

■ regarded as an effective low energy theory valid up to some unknown energy scale and it is an approximation to a more 
] fundamental theory beyond this scale. After having introduced this cutoff parameter e, all the measured quantities 
I that we are considering in the theory and the parameters given in the Hamiltonian become dependent on it. At this 

■ stage, if we remove the cutoff parameter we again encounter the divergent results for the obscrvablcs. However, if we 
, think one of the parameters in the theory (e.g., coupling constant) as a function of e and relate it to an observable 

J> ' (e.g., bound state energy of the system), by solving the appropriate set of equations, we may remove the dependence 
on this unknown scale. That is to say, we can find finite and sensible results for the other observables in the system 
(such as differential cross section, phase shift) by substituting the expression for the coupling constant found in the 
[ previous step and removing e at the end. If all the observables arc still finite after this awkward procedure, the theory 
" " is said to be renormalizable. If not, one must continue to apply the same procedure for other remaining parameters 
(such as charge, mass, etc.) until every observable becomes finite. This renormalization procedure can usually be 
done perturbatively and only few non-perturbative approaches are available since most quantum field theories are not 
exactly solvable. 

When de Broglie wavelength of a particle is much larger than the range of the potential, the interaction can well 
be approximated by a Dirac delta function (point interaction). This problem in one dimension is rather easy and its 
solution is given in any standard textbook in quantum mechanics. If we extend this problem into the one where a 
particle scatters off a periodic set of delta function potentials, it is one of the few completely solvable models Q, which 
describes the electrons moving in a one dimensional crystal lattice. In two and three dimensions, the point interactions 
give rise to infinities but this problem can be cured with the renormalization procedure 0, @, . Most concepts in 
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field theory, such as dimensional transmutation, regularization, renormalization group, etc. can be understood in this 
simpler context. Beside the role that it plays in understanding renormalization, it has many applications in diverse 
areas of physics, as well (see the references in (2I. [lol|). 

Point interactions are also considered in a more rigorous context, so-called self-adjoint extension theory developed 
by Von Neumann and a systematic exposition of this subject has been discussed thoroughly in the monograph [2[, 
where a brief history and an extensive bibliography of it is also given. The formal Hamiltonian in D dimensions 

can be rigorously defined as a self-adjoint extension of a free formally Hermitian Hamiltonian on a space with 
one point removed, where the delta center is located and a boundary condition for the wave function at that point 
is introduced @. Moreover, there is another rigorous approach to the above problem where a relation between the 
resolvents of two different self-adjoint extensions of one symmetric operator is given and it is called Krein's formula. 
The discussion of it for point interactions has been given in [Tl| . Within this formalism, one can immediately 
investigate the spectral properties of the point interactions whereas the domain issues of the operators can be preferably 
handled in the Von Neumann's approach. The results of the self-adjoint extension methods and the renormalization 
approach to the point interactions are the same if a certain relation between the parameter of the extension and the 
renormalized (or bare) coupling constant is satisfied @. 

Many-body version of the point interactions is also extensively discussed in the literature from various directions. 
The Hamiltonian of the system, in which n particles of mass m are interacting through the two-body Dirac delta 
interaction, is 

fc2 " " 

^--£;;Ev?-^ E ^'"^(x.-x,), (2) 



Kj=l 



where A is the coupling constant. One of the earliest studies on the many-body or few-body version of this model 
in two or three dimensions dates back to the work of G. Flamand [l^ and the unpublished thesis of J. Hoppe [l3| . 
and the ones in the Soviet Union, see the references given in 0] . More recently, a perturbative renormalization to the 
above n-body problem has been worked out in [3l and also three-body problem in two dimensions is discussed in 
[Tsj . It has been proved that the perturbative treatment of the three-body problem shows new divergences in three 
dimensions after the renormalization of the two-body sector of the problem and these divergences appear for each 
added new particles Therefore, n — 1 new scales emerge after the renormalization of the rt-body problem. The 
same model is also rigorously studied in [l^ . 

In one dimension, there is no divergence at all and the ground state of this many-body problem is exactly soluble 
p7| and Hartree approximation gives exactly the same results for large values of n [isj . Moreover, the same problem 
for the repulsive case is worked out in [l9| and ^-matrix approach for both the attractive and the repulsive cases has 
been studied in [13, [HI 

A quantum problem where a single particle interacts with a Dirac delta potential in two dimensions shows also an 
elementary example of dimensional transmutation 0, [^, [2H (this term is originally introduced in (23|). Under the 
scaling transformation x — > ax, the Laplacian and i5^^^(a;) function transform similarly. In other words, they have 
the same dimensions [L]~'^ so that the coupling constant A is dimensionlcss in natural units. Therefore, Hamiltonian 
P]) in two dimensions does not contain any intrinsic energy scale due to the dimensionlcss coupling constant. A new 
parameter specifying the bound state energy is introduced after the renormalization procedure, which then fixes the 
energy scale of the system and this phenomenon is called dimensional transmutation. In fact, as shown in Q, the time 
dependent version of this problem has a larger symmetry group 50(2, 1), which exhibits one of the simplest examples 
of anomaly or quantum mechanical symmetry breaking. Furthermore, renormalization group (RG) equations of point 
interactions have been discussed in jl^i H, ESl ■ The /3 function has been calculated exactly in there and the theory 
has been found as asymptotically free in two dimensions. The RG equations for the two dimensional many-body 
extension of the problem, where the Hamiltonian is given by ^ for D = 2, has been addressed for two-body sector 
in (25! . [26j . They are especially useful in this case since there is no analytic solution to the problem. 

S. G. Rajeev [27| introduced a new non-perturbative renormalization method developed for bound state problems 
of some quantum many-body theories: fermionic and bosonic quantum fields interacting with a point source with two 
internal states and non-relativistic bosons interacting via two-body point interactions. One of the main advantages 
of this approach is that all the information about the spectrum of the model is described by an explicit formula 
instead of imposing the boundary conditions on the operator as in the case of self-adjoint extension theory. Another 
advantage is that the renormalization is performed non-perturbatively by introducing fictitious degrees of freedom 
via othofermion algebra so that it helps us to reduce the renormalization to simply normal ordering of an operator 
which is called principal operator $ and then all the information about the spectrum of the problem can be found 
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from its explicit well-defined form. Due to the non-perturbative nature of this method it is also particularly useful for 
dealing with the bound state problems. We are not going to review the original method developed in there. Instead, 
we suggest the reader to read through the relevant parts of the paper [l^], especially ^4'f2+i)NR ™odel to make the 
reading of this paper easier (the problem where bosons interacts with each other via two-body Dirac delta potentials 
is indeed known as the formal non-relativistic limit of the X(f>'^ scalar field theory @, [H, H^). A mathematically more 
rigorous discussion for this approach to X4>(2+i)nr been given in [s^] ■ 

Following the original ideas developed in 27[, we previousl y co nsidered the bound state problem for N point inter- 
actions in two and three dimensional Riemannian manifolds l3l| by using the heat kernel and discussed its spectral 
properties in there. The same model from the Krein's point of view has been discussed for special explicit mani- 
folds, such as strips or tubes [32l. [33j. and it is considered as a natural model for quantum wires including point like 
impurities. The model that we will now construct is the many-body version of our previous work (3l| . where the 
non-relativistic bosons interact with each other via two-body Dirac delta function potential. Our primary motivation 
here is to find a better understanding of the renormalization of many-body models on Riemannian manifolds. 

The paper is organized as follows. In Section II, we construct a model where the non-relativistic bosons interact 
with each other via two-body Dirac delta potential in two dimensional Riemannian manifolds. This construction is 
motivated by the work [27| where a new non-perturbative renormalization method is developed. By extending the 
Fock space, it becomes possible to renormalize the model non-perturbatively by simply normal ordering of an operator, 
called principal operator. Section III is about the mean field approximation of the model and it has been found that 
the magnitude of the ground state energy grows exponentially with the number of bosons, which agrees with answer 
in the flat case already found in The same formulation can also be applied to the one dimensional model where 
there is no renormalization. In this case, the mean field approximation that we develop here gives exactly the same 
result with the one given in the literature p7l [l8| . Finally, we proceed with the renormalization group equations for 
this model and the /3 function is exactly calculated. 



II. CONSTRUCTION OF MODEL 



The non-perturbative method which is applied to our system here consists of the following series of steps: 

1) We first regularize the Hamiltonian via heat kernel 

2) Then we extend the Fock space by the help of orthofcrmion algebra so that new Fock space becomes a direct 
sum of two Hilbert spaces. 

3) The Hamiltonian on the extended Fock space is constructed in such a way that the regularized resolvent projected 
onto the old Fock space gives an equivalent expression for the regularized resolvent of our original Hamiltonian. 
Hence, the coupling constant becomes additive rather than multiplicative. 

4) By normal ordering the equivalent expression of the regularized resolvent, the singular part of the problem 
becomes transparent due to the short time asymptotic expansion of the heat kernel. Then, it is possible to 
choose the coupling constant in such a way that the singular part is removed. 

The Hamiltonian on a two dimensional Riemannian manifold {M, g) is formally given in the second quantized 
language (we use the units such that h = 2m = 1) 




where d^x = y/detg dx^dx^ is the two dimensional volume element, is the Laplace-Beltrami operator (or simply 
Laplacian) defined in a local coordinate system, also written as a; = {x^jx"^) 

and A is a positive coupling constant (it corresponds to an attractive interaction). Here, 0j(.i;), 4>g{x) are the bosonic 

creation-annihilation operators and Sg {x,x') is the Dirac delta function defined on the two dimensional Riemannian 
manifold with metric structure g: 

I dlx5f\x,x')f{x')^f{x). (5) 

JM 
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It is important to notice that the number of bosons d^cc 0j(a;) 4>g{x) is conserved in our model. 

Let us suppose that there exists a negative bound state energy Ei, < corresponding to the normahzed wave 
function tjj{xi, . . . ,Xn',g), that is, 

/ dlxi...dlxn\ipixi,...,Xn;g)f = 1 ■ (6) 

Due to scale invariance of the Hamiltonian under the transformation g — >■ a^g with a positive constant a^, the wave 
fmiction Tp{xi, . . . ,Xn', g) = C("4'{xi, . . . ,Xn',a^g) satisfies the same eigenvalue equation with the energy — a^|£^b|. 
Therefore, the existence of a negative bound state energy implies that the energy can be made arbitrarily negative by 
choosing arbitrarily large values of a. This means that the energy is not bounded from below, which is not allowed 
in a sensible theory. 

In order to cure the problem, we will first regularize the model. The same model in flat space has been discussed in 
[27I Isoj and the renormalization has been performed in a non-perturbative way. In that case, the divergence appears 
due to the large values of momenta (ultraviolet), or short distances. Hence, we expect that the ultraviolet divergence 
must also exist for the same model defined on manifold since every Riemannian manifold can locally be considered as 
a flat space. In [3l|, we have proved that the divergence due to short distance is replaced with the short "time" for a 
simplified version of this model, where a particle interacts with several external delta potentials on a manifold. This 
is accomplished by expressing the resolvent of the system in terms of the heat kernel. In this way, we have been able 
to subtract the divergence from our model by using the short "time" asymptotic behavior of the heat kernel. This 
motivates us that the proper regularization for the many-body version must also be performed via heat kernel and a 
natural choice for the regularized Hamiltonian is 



= Ho-^ f d^xid^x'id^xzd^^d^y <pl{xi)(j>l{x2)K,{xi,y;g)K,{y,x2;g) 

X K,{x[,y]g)K,{y,x2;g)(j)gix[)(t)g{x2) , (7) 

with e the short "time" cutoff parameter, Hq the free Hamiltonian, and K^{x,y;g) the heat kernel on the manifold 
defined as a fundamental solution to the heat equation [s^] 



dKt{x,y;g) _ ^ 



dt 



yiKtix,y;g) . (8) 



Unless otherwise stated, it is always assumed that the Laplacian acts on the functions of the variable x. One of 
the most important properties of the heat kernel that we use in this paper is that it converges to Dirac delta function 

Kt{x,y;g)-^6^^\x,y) , (9) 

as i — > 0+ in the sense of distributions. It is also symmetric Kt{x,y;g) = Kt{y,x;g) for all x,y £ A4 and t > [s^ . 
If we remove the cutoff, that is, take the limit e — >■ 0^, we immediately see that we recover the original Hamiltonian 
given in ([3]). It should also be pointed out that we consider the coupling constant in ([7]) as a function of the cutoff e, 
and its explicit form will be determined later. 

Now, we will consider the resolvent of the Hamiltonian ^ in Fock space J-g with arbitrary number of bosons. 
Following the method developed for the same model in the plane [13] , we will extend the bosonic Fock space J-^ that 
we have started with to a larger Fock space, as it was first introduced in For this purpose, we define new creation 
and annihilation operators, which obey orthofermionic algebra [ssf : 

Xaix)xliy) = <5f (x,y)Ho, 

xgi^)xgiy) = o = xlix)xliy) , (lo) 

where 

Hi=/ dlxxlix)xg{x), no = f-Hi (II) 
Jm 

are the projection operators onto one-orthofermion and no-orthofermion states, respectively. It follows easily that 
there can be at most one orthofermion in any state. The new Fock space is introduced as a direct sum of two Hilbert 
spaces 



(12) 
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where the first sector which docs not include any orthofermion is written as bosonic Fock space and the second 
sector with a single orthofermion as J^g (g)£^(A^). Here, we identify the space of single orthofermion states by £^(A^). 

The advantage of introducing this trick is that it allows us to rewrite the resolvent of the model in such a manner 
that the coupling constant appears additively rather than multiplicatively. Actually, the idea of introducing unphysical 
particles in such a way as to cancel the infinities is not a new idea (see the references in [1^). As a result, we will be 
able to subtract the divergence from our model nonperturbatively by simply normal ordering the operators. Now we 
define the augmented regularized Hamiltonian H'^ on J-^ as 



= HoUo 



V2 



dga:;idga;2dgy (j>l{xi) (l)l{x2)K^{xi,y; g)K^{y,X2; g)Xgiy) + h.i 



Hi 
A(e) 



(13) 



If we split the Hilbert space according to the orthofermion number, the corresponding splitting of the operator 
H"^ — EHq can be written in the following matrix form 



bl 



(14) 



with a : Tb ^ J^B, b\ : Fb ® C'^{M) Fb, J^b ® C^iM) J^b® C^iM). Accordingly, the explicit form of the 
matrix elements of the above matrix is 



bl = 



{Ho - E)Uo , 
1 

71 



dlxidlx2dly <j)l{xi) (j)l{x2)K,{xi,y; g)K,{y, X2]g)Xg{v) 



(15) 



Then, one can construct the augmented regularized resolvent defined as {H^ — EHq) ^ and let us suppose that it is 
of the following matrix form 



R^{E) 



(16) 



Incidentally, the energy E here should be considered as a complex variable. One can find cte, 5^ in terms of a, &g, 
and by a direct computation. This could be done in two apparently different but equivalent ways and the formulas 
were explicitly given in the appendix of [27j . One of the solutions to is 



[a^h\d-H,Y' = ^^ = R%E). 



(17) 



This means that Re{E) projected to J^b is just the resolvent of the operator H^. The other solution for [23| is 



= a ^ + a ^ bl[d^ 



bl] b,a- 



Combining both solutions give 



R'iE) = a, = fl-i + bl [<S>'{E)]~^ b, a'^ 



(18) 
(19) 



where we have defined 



^'{E) = ^- \ [ dlx,dlx2dlydlx'^Alx'^dly' K\{xi,y; g)K,{y, X2; g) 



X K,{x[,y'] g)K^{y' ,x'2; g)xliy) 



(j}g{Xl)(l)g{X2) 



Hn-E 



^Ux[)<j)Ux^2) 



XrAy') 



(20) 



which is called the regularized principal operator, in which the coupling constant is written additively. Now, in 
order to see and separate out the divergent part from ((20)) . we will normal order the operators in ((20)) by using 
the commutation relations of the field operators. For simplicity, we explicitly perform our calculations for compact 
manifolds here, but our result is also valid, in p rinciple, for non-compact manifolds by using a similar method that 
we have done for non-relativistic Lee model |37l . Issj . 
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In analogy with the plane wave mode expansion of the field operators in quantum field theory, one can write the 
eigenfunction expansion of the creation and annihilation operators as 



(21) 



where fi{x]g) is the complete and orthonormal eigenfunction of the Laplace-Beltrami operator [3£ 

-Vg/;(x;5) = (Tifi{x;g) 
dgX fi{x;g)f^{x;g) = Si^n 



(22) 



with the spectrum {0 ~ uq < ai < a2 < ■ ■ ■} so that the free Hamiltonian becomes 



(23) 



It must be emphasized that the degeneracy is formally taken into account in the above sum by the index I. For 
simplicity, we have suppressed this possible degeneracy. Using the commutation relation [<pi,(l)\,] = Sw, it is easy to 
see that {Hq - E)(t)\ = 4)] {Ho - E + ai). Muhiplying this equation by {Hq - E)'^ from left and by (Hq - E + ai)^^ 
from right we get 



(Ho - E)-'<j)l = ^liHo -E + ai)-' . 
We now multiply both sides of the above equation with fi{x; g) and take the sum over / to obtain 



{Ho-E)-'4l{x)^Y.^\ / dte-'(^"-^+'^')/i(x;g) 
, Jo 



(24) 



(25) 



where we have used the fact {Hq - E + ai)'^ = dt e-*(^"-^+'^'). Since (j)} = Jj^ d^y (2/)/;(j/; g) and the 
eigenfunction expansion of the heat kernel is given by [39} 



Ktix,y;g) = Y,<^~'"'M^;9)niy;9) 



(26) 



we find 



{Ho E)-^^\{x) = / dly 4>l{y) / dt e-*(^"-^)/^t(a;,y;5) • 



M 



(27) 



Similarly, by using the same procedure, we can shift all the creation operators 4>\{x'i)4)li{x'2) to the left 

/ d^yid2y^0t(y;)^t(y^) r dte-'("«-^^ Kt{x',,y[-g)Mx',,y'2-.g) 



Hn-E' 



(28) 



and then normal order the new expression with the annihilation operators (j)g{xi)(j)g{x2) so that we obtain the normal 
ordered regularized principal operator 



= ITT " J / 4 (^) / d^^i d2x2 d]x\ di 



X2 </'I(2^l)'/'U2^2) 



X / dt Kt+^{x[,x';g)Kt+e{x',x'2;g)Kt+eixi,x]g)Kt+^{x,X2;g)e ^^gixi) (t)g{x2) 





+ 4/ d2a;id2.T2 </)j,(.Ti) / dtKt+eixi,x';g)Kt+2e{x',x;g)Kt+,{x,X2;g) e-'^"°-^'^q^gix2) 



+ 2 / dtKf+2^{x,x';g)e 
Jo 



~t{Ho-E) 



Xg{x') , 



(29) 
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where the semi-group property of the heat kernel 

Kti+t2ix,y;g) = A^z Kt^{x,z;g)Ku_{z,y]g) (30) 
JM 

is used. We expect that as e — >■ 0+ the last "time" integral in (|29p is divergent since it is the term that corresponds 
to the infinite expression in the principal operator for the flat space R^, where it has been discussed in [27| . In fact, 
we can also naively show that the divergence which appears in the principal operator (|29l) is due to the short "time" 
asymptotic behavior of the heat kernel. 

In order to see this, let us find an upper bound to the expectation value of the last term in the principal operator 
([29)) after taking the limit e 0+. For [n — 2)-bosonic and one-orthofcrmion states 

I d^xxj(x)^(x)|0) , (31) 

JM 

we get for the expectation value 



(*! / dlxAlx' xli^') I AtK^{x,x'-g)e-'^"''-^hs{^m 



oo 



dt (v4"-''|e-*(«°-^)|v4"-'') / dlxdlx' r{x')Kt{x,x'-g) Kt{x,x' ■g)^{x) 



M^ 



<l dt{4-'^\e-'i^o-E)^4-'^) dlxdyKUx,x';g)mxr 



JM^ 

oo 
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<l d>- / dt(4"-^)|e-*(««-^)|^i"-^V^'2t(x,x;g)|V'(x)p 

JM Jo 

poo 

d'^x mxr{4"-''>\ / dt e-'^"°-^^K2t{x,x;g)\4"-''>) , (32) 



where we have used the Cauchy-Schwarz inequality with the semi-group pop and symmetry properties of the heat 
kernel. Therefore, "time" integral in the right hand side of ([5^ is divergent due to the first term in the short "time" 
asymptotic expansion of the diagonal heat kernel, which is given by 

oo 

Kt{x,x]g) ^ -——^^Ukix,x) t"" , (33) 

^ ' k=0 

for any D dimensional Riemannian manifold without boundary j40| . Here Uk{x, x) are scalar polynomials in curvature 
tensor of the manifold and its covariant derivatives at point x G Ai. This means that if the left hand side of (|32|) is 
divergent, this is basically due to the singular behavior of the heat kernel near t = in the last term of the principal 
operator i3S\) . 

All these suggest us to choose the bare coupling constant as 



1 



OO 



A(.) A ''^^ ^''^ 

where — /i^ is to be related to (the experimentally determined) bound state energy of two-boson system. The parameter 
/i^ is at present an arbitrary renormalization scale, which breaks the scale invariance in the unrenormalized problem. 
Even if there is no bound state in the spectrum, our prescription will lead to a finite formulation. Yet, later on we 
will prove that for sufficiently large values of ji^ we can always find a two-body bound state and hence we may solve 
/Lt^ in terms of the physical two-body bound state energy (see equation ([5^ ). In Section |Vl a different prescription 
will be used where the renormalization scale is not directly related to the bound state energy. 

With the present choice of the coupling constant we take the limit e — >■ 0"'" in (P^ . and readily obtain 



<f{E)= j^ dlxdlx'xlix) dt ^_--42)(^^^')_^2(^^^,.^)^-t(/^„-^.)J^^(^^,) 

-\ [ d2xd2x'xj(x) / d^xid^x2d^a;ld^4 0t(a:;)0t(4) / dt Kt{x'^, x' ■ g) 

^ JM^ I J M'^ JO 

X Kt{x' ,X2;g) Kt{xi,x;g) Kt{x,X2;g) e"*^^°"^Vg(a;i) (/>s(a;2) +4 / d^cci dga;2 0j(a;i) 

JM'^ 
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X / dt Kt{xi,x';g)Kt{x',x;g)Kt{x,X2;g)e 



(35) 



This is a well-defined form of the principal operator and we can show that the choice for the coupling constant |34|) 
is sufficient to remove the divergence from our problem. Once we have a proper and well-defined expression of the 
principal operator, we expect that all the divergences are removed since the resolvent which determines the spectrum 
of the problem is expressed in terms of it. It must be emphasized here that the principal operator can be extended 
to its largest domain of definition in the complex energy plane by analytic continuation. 

We must first note that the behavior of the off-diagonal term of the heat kernel near i = is intimately related to 
the small distance behavior due to the initial condition given for the heat kernel. In fact one can show that the choice 
for the coupling constant (j34p is the appropriate one to get rid off the infinity by writing the square of the heat kernel 
in the following subtle way near t = 0: 



dlxd^x'xlix) 



dt 



~87rf 



-S^^\x,x')-K2tix,x';g)Sl;'\x,x')e 
+ Regular terms . 



-t(Ha-E) 



(36) 



The following heuristic argument can be given to justify this choice. Here, what we mean by "regular terms" are the 
other terms in ([55)) and the ignored terms that is coming from the outside of the region t = Q. Let us first look at the 
matrix element of the second term in the first "time" integral in the principal operator (|35p : 



M 



d^a; iPl{x)Kt{x,y;g)Kt{x,y;g)ipb{y) , 



(37) 



as t — > 0+. As a consequence of ([9]), it is possible to replace the function ipl{x) by ipaiv) this limit, so that we have 



d.a; i^aix) Kt{x, y; g) Kt{x, y; g) ipbiv) 



M 



^liv) I <ilx Kt{x,y;g)Kt{x,y;g)iljb{y) 



M 



i^*a{y)K2tiy,y;g)My) > 



(38) 



where we have used the semi-group property of the heat kernel (pO|) . Therefore, if we take the integral (|36|) over x' 
and substitute the first term in the asymptotic expansion p3p of the diagonal heat kernel as t — >■ 0+ , we get 



dgXxlix) 

M Jo 



dt 



-t{Ha-E) 



Stt 



= ^ / d> x\{x) In 



M 



Ho-E 



Xg (x) + Regular terms 



Xg{x) + Regular terms , 



(39) 



where the other terms in the asymptotic expansion p3p do not give rise to an infinite result. 

Let us give a better justification of this choice: we will again assume that the orthofermion operators act on some 
smooth functions; since the set of smooth functions are dense in the Hilbert space norm, this is allowed. We will write 
one of the heat kernels as a distributional solution in ([37]), and use the fact that — is a self adjoint operator. 



d'ydlx raix)Kt{x,y;g)e'^l 5f\x,y)e-<"--^^Uv) 



A12 



M 



e sipl{x)Kt{x,y]g) 



42)(x,y)e-*(^"-^)^,(y) 



Let us expand the exponential e 9 into a formal power series and define 



if fc = 0,2,4,6,.. 



V^(V2)(fe-i)/2^ iffc = 1,3,5,7,..., 



where (Vg/)* = g^-' for any smooth function / on A^. Then we get terms of the following form 

[iV,r~'^Iuix,y;g)] . 



(40) 



(41) 



(42) 
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As t — >■ 0+, the most singular terms in this expansion will come from the terms with the highest number of derivatives 
of the heat kernel, thanks to the following theorem (Lemma 1.7.7 in [13]): If is a differential operator (acting on 

the functions of variable a;) of order a, then the asymptotic expansion of the kernel of the operator D"e*^s on the 
diagonal (in D dimensions) 



fc=0 



(43) 



where Ck are smooth local invariants of the jets of the symbols of the operators and V^. Also arc zero if fc + a 
is odd. Thus, the most singular terms will come from the highest powers of the Laplacian acting on the heat kernel 
when wc formally expand the exponential operator. This means that the dominant contribution to equation (j40p is 
given by 



If wc make use of the heat equation ([8]) in the above, we may infer that 



(44) 



(45) 



Using the fact that e* a*^ generates a time translation by an amount t, which is again true in the sense of distributions: 



lim e* at Kt{x,x';g) = lim Kt+t'ix,x'; g) = K2tix,x';g) , 

t'—yt t'^-t 

we see that the most singular part of the integral as i — > 0+ turns out to be 

dly raiv) K2tiy, y; g) e-*(^°-^) Mv) , 

M 

where we have taken the integral with respect to x. This justifies our choice of the coupling constant ([34]). 



(46) 



(47) 



A. Two Dimensional Flat Case Revisited 



Wc can also explicitly show that this idea works for the same model on flat space M? by writing the principal 
operator in momentum space that has already been calculated in 27[. For this purpose, let us consider the first part 
of equation p5|) in a two dimensional plane, i.e.. 



d^xd^a;'x'''(x) lim^ / dt 



-tfj,'' 



■(5(2)(x,x')-i^t^(x,x')e 



/^„-^(//o--B) 



X(x') . 



Substituting the explicit form of the heat kernel in [34 



ift(x,x') 



-|x-x'|^/4t 

Ant 



we find for El 



(48) 



(49) 



d^xd^x' xH'^) lim / dt 



-tfJ.^ 



SlTt 



-,5(2) (x,x') 



Hx-x'[V4t g-|x-x'lV4t 

(Ant) (Ant) 



X(x') . (50) 



.2 

If we write the heat kernel as a Fourier transform of a function p and then change the integration order above for 
the second term, we get 



d^xd^x'xH^) / 



= -|x-x'|^/4t g-|x-x'p/4t 

(47rt) (47rt) 



-t{Ho-E] 



X(x') 
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dt 



gip.(x-x')-tpV2 

(2^ (8^rt) ' 



-t{Ho-E) 



X(x') 



(2' 



dt 



i&TTt) 



-X(P) 



Then, equation p8)) becomes as e — > 0+ 



dt 



87rt 



^-t{Ha-E+p^/2) 



go-g + pV2 



This is exactly the same result that was already calculated for this model defined in the flat space . 



(51) 



X(P) • (52) 



B. Analysis of the Bound State Spectrum 



As a result of our analysis, we now have obtained a finite well defined model, that is, the resolvent of the system is 
expressed in terms of the well-defined principal operator given in ([35]) 



R{E) 



1 



1 



1 



Hn — E 2 Hn 



E 



dlv <i>\{y)^\{y)xM <^~He) / d^y Uy)Uv)x\{y) 



1 



Hn-E 



(53) 



IM J M 

This is the analogue of the Krein's formula in the case of the many-body version of the point interactions. All the 
information of the spectrum of the problem can be determined from the above resolvent operator. In this subsection 
we will discuss the spectral properties of our model, especially bound state spectrum. 

The poles in the resolvent corresponds to the bound states. For non-compact manifolds, there can not be any pole 
due to the free resolvent. For compact manifolds, we are interested in the poles below the poles of the free resolvent. 
These imply that the roots of the principal operator ([35]) 

^{E)\^) = , (54) 

determines the possible bound state spectrum. As in the case of the problem where the particles onl y in teract with 
an external Dirac delta potential, which displays a dimensional transmutation in two dimensions [3l[5ll22l|. our model 
constructed above also realizes a kind of dimensional transmutation. This can be seen as follows. From the original 
Hamiltonian ([3]) that we have started, it is easy to see that the coupling constant is dimensionless so that there 
seems to be no parameter whatsoever to yield an estimate of the energy by naive dimensional analysis. However, 
if we have a length scale coming from the geometry, such as the curvature, this provides a geometric energy scale 
which is there also for the free theory. Nevertheless, even if it is the case, a new dimensional parameter /x^ shows up 
after the renormalization procedure from the relation (j34[) . Therefore, we can say that this is a general dimensional 
transmutation and it is most striking when there is no intrinsic energy scale coming from the geometry [3]| . 

After the renormalization of the coupling constant, we must be able to predict the other measurable quantities in 
terms of the measured two-particle bound state energy E^^J , in our version the arbitrary scale — /i^ should be solved 
in favor of this binding energy. In flat space M^, two-body solution is given by iSgr"* = — /^^ WA - From this point on 
we assume —fp is expressed in terms of Ejj^K We make the following comment, let us consider a compact manifold, 
and apply the variational principle for the first eigenvalue ujq{E) of $(i?) in the two-boson sector. Since we are on a 



compact manifold we choose the orthofermion wave function as constant, 
value of the principal operator $(i?) by the following variational ansatz 



= |0) ® 



7m) L"'-''"''"^^ 



We now calculate the expectation 



(55) 



Since $(£') is normal ordered, all the parts which contain bosonic creation and annihilation operators will vanish. 
The only term which survives sets an upper bound for ujq{E). Hence, 



< / dt 
Jo 

/•OO 

< dt 
Jo 



Sirt \V{M)Jm^ ' ' 
e.-^'* / 1 



dlx dlx' Kt{x, x'- g)Kt{x', x; g) ] e"!^!* 



\V{M)Jm ' 



d X K2tix,x;g) 



-\E\t 



(56) 
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where we have used the semi- group property of the heat kernel (|30)) . Compactness of the manifold implies that it is 
complete as a Riemannian manifold and it has a Ricci tensor bounded from below which we formally write Rc > k. 
As a result of the theorem proven by J. Cheeger and S.-T. Yau [4l|, the heat kernel has the following lower bound 

Ktix,y;g)>Kt^idg{x,y)), (57) 

where is the heat kernel of the simply connected complete two dimensional manifold of constant sectional curvature 
K. In particular, we choose K^{dg{x,y)) as the heat kernel of the two dimensional Hyperbolic manifold H"" for 
K = —1/R?, where R is the corresponding length scale. In case the lower bound is positive we may choose the heat 
kernel for the two dimensional flat space and the argument below becomes even simpler. Since the heat kernel for 
two dimensional Hyperbolic manifolds is explicitly known j42|, a lower bound of the diagonal heat kernel in ([55)) is 



From the expansion of the function cosh, we can write the denominator as V cosh s — 1 = 5^*^/(2^:)! 
is/V2)y^J2'Zi 2.s2'^-V(2fc)!- Then we have 



t^o{E) < dt 







,(59) 



where we have added and subtracted 1 in the parenthesis above. Since 

2s2fe-2/(2/fc)! < for ah s > 0, we have 
u;,iE) < ^ In f ^^-tl) + Tdt-Lfi- e-*m+^) . (60) 




Using VT+4t < 1 + 2i and 1 + 2t > 1 for all t > 0, we get 

u^oiE) < ^ In f ^^-^ I + ^-—^ ■ (61) 



(62) 



For large values of fj,^ there always exists a unique E^, < such that 




; 2tt\E,\ + ^ 



As we will prove in this section 



dujo 
dE 



< , (63) 



to get the true zero Egr of ujq{E), we must further decrease E (or increase \E\) so that we will have a well-defined 
expression of /i^ in terms of two-particle binding energy Egr^ < < 0, as shown in Figure [T] Therefore, by assuming 
that two-body problem is solved, we can then study n-body problem. Since we are only interested in the bound states of 
the model at the moment, we should be able to determine 7i-particle bound states after the renormalization procedure. 
The exact treatment of this problem is rather difficult. Assuming that the details of the two-body interaction can be 
understood, we will study the model in the mean field approximation in Section IIIII Before embarking on studying 
the mean field analysis, we will make some general remarks about the bound state spectrum of the problem. 

It is a well known fact that the residue of the resolvent at its isolated pole /i is the projection operator to the 
corresponding eigenspace of the Hamiltonian 



— (f dE R{E) , (64) 



27^^ 

where F^ is a small contour enclosing the isolated eigenvalue /x in the complex energy plane [43|. Let us suppose 
that there exist a ground state and choose our contour enclosing this ground state energy, namely Eg^. Then, the 
above integral of R{E) gives the projection to the eigenspace |^o)(^o| corresponding to the minimum eigenvalue of 
the renormalized Hamiltonian for the many-body system. 
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upper hound 




FIG. 1: A Typical Flow of the First Eigenvalue of Principal Operator. 

From ([35]). it is easily seen that the principal operator formally satisfies = ^{E*). We assume that ^{E) 

defines a self-adjoint holomorphic family of type A [i^l, so that we can apply the spectral theorem for the principal 
operator or inverse of it. Since the principal operator <^{E) acts on J^g^""^' ® C'^{M), we have 



^-\E)=Y, 



¥u{E) + / ME) 



u:{E) 



(65) 



where the projection operator 



n{E) = \<Pk{E)){(j)k{E)\ = \LOk{E)-ME)){ujk{E)-nk{E)\ , 
is given in terms of n — 2 bosonic particle state and one-particle orthofermion state: 



(66) 



dixMx)xiixm ■ 

M 



■X„-2) 



(67) 



Here uJk{E) and \LL!k{E);flk{E)) are the eigenvalues and the eigenvectors of the principal operator, respectively. 
Similarly, the (generalized) projection operator 



F^E) = \<f>iE)){HE)\^HE);n{E)){L,iE);niE)\ 



(68) 



corresponds to the continuous eigenvalues and eigenvectors of the principal operator. We assume that the principal 
operator has discrete as well as continuous eigenvalues and the bottom of the spectrum corresponds to a non-degenerate 
eigenvalue. Above integral is taken over the continuous spectrum cr($) of the principal operator (for simplicity, we 
write it formally, it should be written more precisely as a Riemann-Stieltjies integral). 

As emphasized in the previous section, the bound state spectrum corresponds to the solutions of the zero eigenvalues 
of the principal operator (j35p . In order to estimate the ground state energy of our system, it is crucial to determine 
how the eigenvalues Wfc evolve with E. For this purpose, let us calculate the derivative of the eigenvalue ujk of the 
principal operator with respect to E. If we apply the Feynman-Hellman theorem to the eigenvalue problem for the 
principal operator, we get 



dE 



, dHE) 

' dE ' 



(69) 



A direct computation for the derivative of the principal operator ((35|) with respect to the energy E gives 
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dlxdlx'xlix) r dtt K^{x,x';g)e-'^""-''\g{x') + ^ [ d^x d^x' d^xi d^X2 d^x'i dg4 



X xlix) Xg{x') 4'l{x\)<t>l{x'2) / dt t Kt(xi,x; g) Kt{x2,x; g) Kt{x' , x'^-g) Kt{x' , x'^- g) e-*(^°-^)(/.g(a;i) ^^(xa) 

Jo 

+2 / d^xd^x dlxidlx2xl{x)Xg{x)4'l{xi) 



X / dttKtix2iX;g)Ktix,x';g)Ktix',Xi;g)e-"^"°-'^Ug{x2) 



(70) 



For simplicity, we will separate the terms in the expectation value of the principal operator in (|69p . using (j70p . Let 
us first consider the first term 



d^xdyrix) / dttKf(x,x';g){ujk\e~'^"°''^^\iJkH{x') 
A42 Jo 



(71) 



where ip{x) is the wave function of the orthofermion. If we think of the factor t in the above integrand as an integral 
J_j(du/2) and then make the change of variables t = ti + t2, u = ti — t2, we readily obtain 



00 /'OO 



dlxdlx'rix) / / dt,dt2 ir,%,,(x,x';g)(wfe|e-(*^+*^)(^°-^)|^fc)V(x'). 
Jo Ja 

Using the semi-group property of the heat kernel pop , equation (|72p can be rewritten as 



(72) 



M- 



/•OO /"OO p 

d^xd^gx' tp*{x) / dtidt2 / dlzidlz2 Kti{x,zi;g)Kt^{zi;x';g)Kti{x,Z2;g) 

Jo Jo J 



X A'4,(z2,x';g)(c.fc|e-(*i+*^)(««-^)|L.fc)V'(a;') 



(73) 



Changing the order of integrations we find 



dgZl dgZ2 



g-ti(ffo-iJ) / dlxKt,{zux;g)Kt,{x,z2;gW{x)\u:k) 



M 



(74) 



which is obviously always positive. Now we return to the expectation value of the second and third terms in (|70p . It 
is easy to see that they can be expressed as 



dt t 



dlxi dlx2 d^x Kt{xi, X] g)Kt{x, X2; g)ip* (x) 



liHo-E] 



4>g{xi)(j)g {X2)\(^k) 



2 dp I dtt 

IM Jo 



dlxdlx2 Ktf2{z,x; g)Ktix,X2; gWix) e-^^"°-^Ugix2)\L0k) 



(75) 



where we have used the fact that we can rewrite the second heat kernel Kt{x, x' ; g) in the third term of (j70p as 
Jj^ dgZ Kt/2{x, z; g)Kt/2{z,x'; g) by the semi-group property (pO)) . Consequently, we obtain 



dujk 
'dE 



< . 



(76) 



The eigenvalues ajfc(£')'s flow with E in accordance with (|76p . that is, these are monotonically decreasing functions 
of E. For sufficiently small values of i?, there can not be a zero eigenvalue of the principal operator since the energy 
must be bounded from below. Moreover, for a given the eigenvalues can be ordered as (-£'*) < wi(i?») < . . .. 
Therefore, due to (j76p and non-degeneracy of the lowest eigenvalue ljq, only the minimum eigenvalue luq flows to 
its zero value at the minimum energy E = Egr. Hence, the ground state corresponds to the zero of the minimum 
eigenvalue ujo{E) of ^{E), as shown in Figure [2j 
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FIG. 2: A Typical Flow of the Eigenvalues of Principal Operator. 

We may now show that E'g"'' < Eg^'^ for compact manifolds. To see this, we take the solution of two-body ground 
state Ivf^^)) as 



^ |0) ® / dgx 4'\x)xHxm 
Jm 

and then make a new ansatz |vl/"'^'"^ for the n-body problem in the form 

It is convenient to split the principal operator given in ([35]) as K{E) — U{E), where 



K{E) 



AlxAlx'xl^x) 



dt 



-6^^\x,x')-KUx,x';g)e~'^"° 



-E) 



(77) 



(78) 



(79) 



and 



U{E) = \[ At ( dlxdlxidlx2Kt{xi,x]g)Kt{x,X2\g)(j)g{xi)(t)g{x2)Xg{x) 



-t(Ho-E) 



dlx' dlx[ dlx'2 Kt{x[,x'; g) Kt{x' , x'2; g) (t)g{x[) 4>g{x'2) Xg{x') 



2 1 d^z 

M 



dt 



dgX d^xi Kt{xi,x] g) Kt/2{x, z; g) 4>g{xi)xg{x) 



d^x' dlx[ Kt{x[,x'; g) Kt/2{x' , z; g) (j)g{x[)xg{x') 



(80) 



where we have used the semi-group property of heat kernel pop and the assumption that we can interchange the order 
of integrations. By the variational principle, 



< (*™n$(£;^?.^)i*™'') = -(*™n(7(4?)i* 



(81) 



where {ip'^'^\K{Egr')\'il}]^') = 0. In order to calculate the above expectation value, we first show that 



n(2)MJ,(2) 



e 2^ 



/ dlx' d^xl dlx'^ Kt{x[,x';g) Kt{x' , 4;^) M^'i) Xs(^')l*™^) 



^ (n-2)V2(n-3)V^ _.|^(., 

[l/(Al)]("-2)/2 



e -"-I / dlx'dlx[dlx'2Kt{x[,x';g)Ktix',x'2;g)4'\x') 
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M 



[y(7W)]("-2)/2 ^ /. .S^ ^0 (2; j 



dgys • • ■dgy„_2|y3, • • ■ ,yn-2) 



(82) 



A4" 



where we have used the fact that the free Hamiltonian operates on bosonic states and gives zero for constant wave 
functions. We have also used the stochastic completeness of the heat kernel in the last line. Then, the expectation 
value of the first term in ([50)1 becomes 



(n-2)(ri-3) 



At e-*l4^ 



M 



[V{M) 



(n-4) 



which is finite due to 



M 



< 



M 



M 



dlx 



[y(7W)]("-2) 



V{M) 



Similarly, we can calculate the expectation value of the second term in (j80p and get 

2' 



2(n- 2) 









/ dt 




j 


'o 


JM 


J M 



dixKt/2{x,z;g)'iJjP{x 



(83) 



(84) 



(85) 



It can be shown that the above integral is finite if we use the eigenfunction expansions ([22)) and pG]) , so that we have 



M 



M 



dlxKt/2{x,z;g)% ' (x) 



dlz 



M 



(86) 



where 'ip^\x) = J2i fiix, g)^(^\l)- Since J2i e"**^' |V'o^^(OP — J2i I'^o^'lOP ^-i^d the minimum eigenvalue <7o = 
for compact manifolds, the above integral is bounded from above by 1 and so that equation (j85p is finite. Hence, we 
obtain 



(n),p(2)N < (n~2)(n~3) 



M 

2(n- 2) 
V{M) 



poo 






/ dt 






'0 


JM 


J M 



d^^x Kt/2{x, z; g)ipP (x) 



e-*K'l<0. (87) 



As a consequence of (ffS)) and (|57|) . to find the zero of ujq{E) in the ri-particle sector we must reduce E below Ejp.\ as 
shown in Figure [3] This completes the proof. We will now calculate n-particlc ground state wave function in terms 




FIG. 3: A Typical Flow of the Eigenvalues of Principal Operator in the two and n-boson sector. 



of the solution |(^o) of ^iEgr)\4>o) = 0- 

Let us expand the minimum eigenvalue u)q{E) near the bound state energy Eg,. 
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UJo{E) = UJo{Egr) + {E- Egr) 



Using this result and (|65l). equation ([64)) yields 



dE 



{E - Egr) 



duJo{E) 



dE 



i(Fo - Egr)-' dlx cpl{x)4>l{x)M-) ( 



j \uj^{Egr)){u:,{Egr)\ 

x[ dlycPg{y)cPgiy)ro{y){Ho-Egr)-' . (89) 
Jm 

We assume that there is no other pole coming from {Hq — E)~' near Egr, and no other terms for fc ^ contribute to 
the integral around E = Egr- Let the eigenvector of the principal operator corresponding to the ground state be 



|0o(-Egr)> = / dga;i • • •dga;„_2 uo(.Ti, • • • ,a;„_2)|a;i • • •Xn-2) / dg.T ?/>o(x)xJ(2;)|0) 



(90) 



M 



By using the eigcnfunction expansion of the creation and the annihilation operators and their commutation relations, 
we will shift all creation operators 4>l{x) in (|M|) coming from ([50)) to the leftmost 



-.'l^lix)(l}lix)(l)lixi) ■ ■ ■ 0],(a;„_2) 



dgyi---d^2/n (jyUyi) ■ ■ ■ (jyUyn) 



Ho — E - " - - j_f^„ 

/>oo 

X / dt e-'^"^-''^ Kt{yi,x;g)Kt{y2,x;g)Kt{y3,xi;g) - ■ ■ Kt{yn,Xn-2;9) , (91) 



and all annihilation operators (f>g{x) in (|5^ coming from (|M|) to the rightmost 



't'g{x)4>g{x)4>g{xi) ■ ■ ■ 0g(a;„_2) 



rr p / dyi---dyn 



dt 



X Kt{yi,x; g)Ktiy2,x] g)Kt{y3,xi] g) ■ ■ ■ Kt{yn,Xn-2; g)(t>giyi) ■ ■ ■ (l)giyn) , (92) 

which are the generalized versions of equations we first used in f37| . Therefore, from equation ()89p . we read the state 
vector 1^0) of our many-body system in terms of the eigenstate \(j)o) of the principal operator 



l*o> 



da2/l • • -dgi/n *o(2/l, ■ • ■ ,2/«)|yi ■ ■ • 2/«) 



/ dgXi • • • dgX„_2dg2; ^ V / dt e *l^'"-lis:t(2/^(i),x;5)iirt(y<^(2),a;;g) 
Ja^"-i Jo 

9wo(^) 



<Te[i---Ti 

X Kt{y„(3,),xi\g) ■ ■ ■ Kt{y„(n),Xn-2\g) uq{xi, - ■ ■ ,Xn-2)ipo{x) 



dE 



-1/2 



\yi---yn), (93) 



where the sum runs over all permutations a of [123 . . .n]. Comparing equation (j90p and equation (|93p. we see that 
the state j^'o) is a complicated convoluted integral of the eigenstate |0o) with the heat kernels. 



III. MEAN FIELD APPROXIMATION 



In standard quantum field theory, one expects that all the bosons have the same wave function u{x) for the limit of 
large number of bosons, i.e., as n — > 00 and the wave function of the system has the product form of the one-particle 
wave functions. However, due to the singular structure of our problem, the wave function in (j93p can not have a 
product form in the large n limit. In order to see this, we scale t = t' /\Egr\. With a hindsight coming from the proof 
that the lower bound of the ground state energy grows exponentially with the number of bosons in flat space [27l . Isoj 
we may assume that Egr grows fast enough as n increases. In this case, all integrals of the heat kernels are peaked 
around ya{k)- (This is clear from © and also from the stochastic completeness assumption). Then, all integrals of x; 
are 
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M 



d^xi Kt/\E^,\ixi,ya{i+i))uo{xi, ...,xi,.. . ,a;„_2) ~ uo{xi, . . . ,ycr(z+i), . . • ,a;„_2) , 



(94) 



for ^ = 1, . . . , n — 2 as n — > oo and similarly for x integral. Then, the state |*I'o) becomes 



If 1 r°° 



't\Ea 



cr6[l---n] 



X Wo(2/<t(3),--- ,y<T(n))V'0(ya(2)) ( 



-1/2 



|2/1 • • • Vn) ■ (95) 



One can understand the singular nature of the wave function in this limiting form more easily. We pick any two 
bosons, and transform them through our formalism into an orthofermion, with its wave function tpQ to be determined 
consistently. This orthofermion wave function corresponding to the pairing, could be quite regular, yet its multipli- 
cation with the heat kernel, integrated over the time variable produces a function singular as the two variables of the 
heat kernel approach to one another. This singularity is the same as the singularity of the bound state wave function 
of a particle interacting with a delta source [3l[, hence it is square integrable. 

It is important to notice that \'^o) is not in the domain of H^. To prove this, it is sufficient to consider the following 
term which appears in calculating (^'o|i?o|^o) 



dlxdly / dti e-*il^-lx,^(a;,y;5)Vo(y) 

/•OO 

dlxdly / diie-*il^-l/C*,(x,y;5)^o(y) 



dt2 e-*^l^« 



1 



'^lKt^{x,y;g) 



dt 



2 e 



2m, 

-t2\E,r\{ _ 9Kt2ix,y;g) 
dt2 



(96) 



where we have used the fact that the heat kernel satisfies the heat equation ((8)). After applying the integration by 
parts to the t2 integral and using the initial condition for the heat kernel Kt{x,y;g) — 5g{x,y) as i — )■ 0+ and (|30p . 
we find 



d^xd^y / dti e-*il^-lA't,(x,2/;g)V'o(2/) 



/•OO 

5g{x,y)-\Egr\ / dt2 e-*^l^-lA'i,(a;,j/;5) 
Jo 



M 



d'x / dtie-'^\^o-\Kt,{x,x-g)^o{x) 







\E, 



gr\ 



dti e 



— t\ \Eg 



dlx 



'0 JM 

After the change of variables u = + ^2 and v = ti — t2, we get 



dt2 e-*^l^-l/^t,+t, (x, x; g)Mx) ■ (97) 



n poo nOO 

/ d^xMx) / dtie-''\^^''^Kt,ix,x;g)-\Egr\ dw u e-"l^-l/C„(a;, x; g) 

JM Jo Jo 



(98) 



The first term is divergent due to ([33]). Similar to the problem with point interactions on manifolds which we studied 
in (3l| . our problem here can also be considered as a kind of self-adjoint extension since the state ^'o does not belong 
to the domain of the free Hamiltonian. The self-adjoint extension of the free Hamiltonian extends this domain such 
that the state ^o is included. Although the state ^'o is not in the domain of Hq, the eigenvector corresponding to 
the eigenfunction uo{x) for the lowest eigenvalue of $(i5) can be taken in the domain of Hq. 
As a result, I'I'o) given in (|95)) is not in the product form in the large n limit, that is. 



dlyi---dlyn Y['fo{yk)\yi---yn 

fc=i 



(99) 



The solution takes a kind of convolution of the wave functions in the domain of Hq with the bound state wave function 
which is outside of this domain. 
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Yet, $(_B)'s lowest eigenfunction may well be approximated by a product form for large number of bosons, that is, 

Uq{xi,--- ,Xn-2) = Uoixi) ■ ■ ■Uo{Xn-2) , (100) 

with the normalization 



\\uo\ 



d'x \uoix)\' = 1 , / d'x \Mx)\ 



1 . 



I M J M 

Therefore, the expectation value of the principal operator by applying the mean field ansatz must vanish, that is. 



(101) 



(0o|$(^g.)l'/'o) = 0. 



(102) 



Although such a mean field approximation is expected to be crude in less than three dimensions, F. Calogero and 
A. Degasperis [l8| have shown that even in one dimension the mean field approach to this problem gives an excellent 
agreement with the exact result. This is a finite problem and we will see in the next subsection that the present 
approach is also consistent with the exact result. Therefore, we expect that the mean field approximation to this 
problem in two-dimensions is also reliable. 

In order to calculate (jl02p explicitly, we will make normal ordering of the creation and the annihilation operators 
by using their eigenfunction expansion. Hence, the equation above yields 



dt 



Snt 



dlxdly 'ip^{x)'ipoiy)K^ix,y;g)e- 



dVdV ul{x')Kt{x',y'-g)u^[y') 



ri-2 



(n-2)(n-3) 



dt 



^xd?xi ^X2 u*^(x{)u*^{x2)Kt(x-Y,x;g) Kt[x2,x;g)'\\}^{x) 



X e f / dlx' dly' ul{x')K,{x', y'; g)u^{y') 



ri— 4 



+ 2(n- 2) / dt 
'o 



dlxdly ul(x)Kt{x,y]g)i)o{y) 



-t\E,A 



XI / dlx dly' ul{x')Kt{x' ,y']g)u(i{y') 



n — 3 



(103) 



We are going to approximately solve Egr from the above equality for large values of n. In order to solve it, we may 
assume that \Egr\ grows rapidly with n. This is plausible because \Egr\ ~ [i^e^^^l^ for fiat space given in the 
mean field approximation j27| . Every Riemannian manifold can locally be considered as a fiat space, and the infinity 
appears due to the high values of momenta (ultraviolet divergence) or short distances we expect that the result for 
the large n behavior of the ground state energy is similar on the manifold case. This allows us to consider the above 
equality in the large values \Egr \ ^ so our aim is to find only the terms that contribute most to the above integrals. 
We first calculate asymptotically the left hand side of (|103p 



dt 



Snt 



dlxdly 'ip^{x)ipoiy)Kt{x,y]g)e 



dlx' dly' ul{x')Kt{x',y'-g)m{y') 



(104) 



for the large values of \Egr\. We will now ignore the additive constants to n, e.g., n — 2 ~ n since n is very large. The 
major contribution to the above integral for large values of \Egr \ can be computed since the asymptotic expansion of 
the following form, namely Laplace integrals 



dt f{t)e 



\E,r\g{f) 



(105) 



is given by Watson's Lemma |45| . The main contribution to the above integral can be obtained by Taylor or when 
necessary by the asymptotic expansions of the functions f{t) and g{t) near the minimum of g{t). Similar to the 
reasoning given in the previous section, we write the square of the heat kernel in a subtle way, that is, we will use the 
initial condition for one of the heat kernels near t = Q. After this and an integration, we substitute the asymptotic 
expansion psp for the diagonal heat kernel near t ~ Q (the region that gives the dominant contribution). Hence, the 
left hand side for large values of \Egr \ ^ /i^ becomes 
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dt 



SlTt 



M 



dlx |Vo(x)| = 



M 



d^x' \uq{x')\' 



dt 



87ri 



Sirt 



^—h-,{\E,r\/^i^ 
on 



(106) 



As for the right hand side of (|103p , we apply the same method while we keep the next order terms coming from the 
eigenfunction expansion of the heat kernel in the n-th power of the integrals. Therefore, we obtain 



2 /-txj 



dt 



dlx |uo(2;)PV'o(a;) 



M 



2n dt 
Jo 



dlx UQ{x)')jjo{x) 



M 



dlx \uo{x)\^ - tK[uo]] , (107) 



M 



where we have defined 



K[uo] = I dlx \VgUQ{x)[' 

J M 



(108) 



and used the eigenfunction expansion of the heat kernel (|26p and expanded the exponential inside by keeping the first 
two terms: 



Kt{x,y;g) tcri) fi{x; g)fi{y; g) 



(109) 



We can rewrite the above expression (|107p by making a change of variable t = t' /\Egr \ as 



2 /-oo 



dt' 

2 Jo \Egr\ 



dlx |uo(a^)|Vo(a;) 



M 



^1 \ l-Egr 

K[u„] 



\Enr 



2n 



dt' 

\Egr\ 



dlx u*q{x)-)Po{x) 



M 



^1 \ \Egr\ 



Moreover, we can think of terms in the square brackets as an exponential for large values of \Egr\ so that 

2 



2 /-oo 



dt' 



n 

2 Jo \Egr\ 



dlx \uo{x)fil;Q{x] 



M 



+ 2?! 



dt' 

\Egr\ 



dqX ul{x)'4)o{,x] 



M 



(110) 



(111) 



From equation (jl03p . it is easy to see that the left hand side is a monotonically increasing function and the right hand 
side is a monotonically decreasing function of \Egr \ so there is a unique solution, say at Below this point ji^grl, 

the left hand side is always less than the right hand side. Therefore, if we can find an upper bound to the right hand 
side of (|llip . and find a solution at | this implies that Egr > — For this reason, let us first set the normalized 
wave function of the orthofcrmion to saturate the Cauchy-Schwarz inequality (as noted similarly in the flat case j27| ) 



■00 (a;) 



1/2 



(112) 



Then, the upper bound of the right hand side of (lllip is 
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1 /■ 9r) 1 I •'■M S 



^ 2 
2 ' 



2n 1 



< TTT^ — r — r- / d„a: Iwofa;)! + tt^ — r t^t — r- , (113) 



where the Cauchy-Schwarz inequality in the second term is used, that is, 



XMd^2;u5(a;)|uo(a;)p) (^|||uo(a;)P|| ||Mo(a; 



< ' r .^M4 ' = 1 ■ (114) 



Im l"o(a;)|4 d^x |Mo(a;)| 

We now recall the following theorem (Theorem 2.21 in |4^): The Sobolev imbedding theorem holds for a D 
dimensional complete Riemannian manifold M with bounded curvature and injectivity radius i5 > 0. Moreover, for 
any e > 0, there exists a constant Aq{e) such that every Lp £ H'l{M.) (HKAi) is the Sobolev space defined on a 
manifold M) satisfies 

Mp < iKiD,q) + e)\\Vgip\\g + Ag{e)M, , (115) 

where 1/p — 1/? — 1/-D and 



D-q\D{q-l)J \T{D/q)T{D+l- D/q)LOD- 

with ujd-1 is the volume of of unit radius. 

Furthermore, there is an optimal inequality for the two dimensional case given by T. Aubin (46l . l47j and it states 
that: Let be a D dimensional C°° Riemannian manifolds with injectivity radius (5 > 0. If the curvature is constant 
or if the dimension is two and the curvature is bounded, then Aq{0) exists and every if E Hl{M) satisfies 

||^||p<K(i^,q)||V<,^||, + A,(0)||^||, . (117) 

For and H^, the inequality holds with A,(0) = 0. 

Let us choose p = 2, q = 1 and D = 2 for our purposes, the inequality (jll7p is reduced to 



1/2 

dlx\^{x)A <- f dlx\VM^)\+Af d^2;|^(x)|, (118) 

M / JM JM 

where K(2, 1) = 2/t: and Ai{0) = A. If we set ip{x) = \uo{x)\'^, then 

dlx\uo{x)A <-[ dlx\u*{x)VgUo{x)\ + - [ dlx\ua{x)VgU*o{x)\ + A [ d^a; |Mo(a;)|' 

M J ^ JM JM JM 

< A + {4/7r)K^/''[uo] , (119) 
where we have used Cauchy-Schwarz inequality and the normalization of uq{x). Hence we obtain an upper bound for 



A+{i/TT)K^/^[uo]] 



2\Egr\ (1 + ^) \Egr\(^l + nigB>\ 



2n 1 , , 

(120) 



-Eg,, I ' ^ l-Egr 

Finally, combining the two results, we find that 

( ^ + FS7.^ 

\Enr 



1 + I3z 

47r "-^'-y'"- ' l + az 



^'^'ln(|£;,,|/M2)<^2^2A__Z_^ (121) 
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where a = \l\Egr\; P = ^/{T^A^/n) and z = ■\/nA'[wo]. For simplicity we ignore the second term in the right hand 
side but we will return to these issues once we find the solution and check the consistency of the approximations that 
we have made so far. An upper bound of the right hand side is achieved at = /3/a and its value is 'n?A^{l + ^). 
As a result of these, we eventually obtain 

Egr > -M^e"(2V^) . (122) 

We note that the location of this maximum for the variable z is only formal, and does not correspond to the physical 
value of nK[u{)\. It is simply chosen to get an upper bound for the right hand side, thus a lower bound for the energy. 
In fact, to be physically consistent, nK[ui3\ should be of the order of \Egr\ in the mean field approximation. Since 
we do not know a method to solve these equations, it is not possible to calculate the actual values. Yet it is easy 
to check that in the limit where nK^ao] 3> |£'gr|, the renormalized term becomes dominated by this kinetic term, 
and the potential part also becomes much less than the renormalized term, hence there cannot be a zero for the 
operator ^{E) under these assumptions. Hence, we can keep iir[wo] ^ \Egr\ condition in our approach. This has a 
nice interpretation physically, for the $(£') operator, the ordinary total kinetic energy is of the order of the binding 
energy, moreover, the binding pair, transformed into orthofermion, has also finite kinetic energy. Nevertheless, we 
know that the actual wave function has infinite kinetic energy, thus this formalism nicely takes out these pairs and 
converts them into regularly interacting particles. As a result, they satisfy a nonlinear eigenvalue equation. 

After wc find the solution, it is easy to check the approximations that we have made, the order of all these ignored 
terms are indeed small. To be more precise, the next order terms coming from the asymptotic expansion become 
lower order terms in n for the ground state energy. 



IV. CONFIRMATION OF THE PRESENT METHOD IN ONE DIMENSION 



We can apply our method to the ground state for the same system in one dimension, where there is no need for 
renormalization as can be easily seen from the short "time" asymptotic expansion of the heat kernel ([55)) in ([5^ . The 
exact solution and the Hartree approximation (for bosons) to the ground state in one dimension have been studied in 
[13, [il). The exact solution is given by [Ol 

^'(a;i,...,a;„) = Cexp j ^ , (123) 

\ i>J=l 

where the normalization condition (Jjg„ dxi . . . d.T„(5(a;c.rn)|^'P = n) allows us to calculate the constant C explicitly 
[l7{ . The exact ground state energy is then 

Egr^-^nin^-l). (124) 

The Hartree solution to the ground state wave function (except for the infinite degeneracy due to translational 
invariance) of the same system [Tsl is 

. , (125, 

cosh [XnxiA) 

where dx ["(^(a;)!^ = 1. Since n is large in this approximation, we may also write the above solution as 
(An/2)^/^e~'^"l^l/* and the ground state energy is 

Ef, = -^n\n-l). (126) 

It is obvious that the exact results for the ground state coincides with the results given in the Hartree approximation 
in the large particle number limit. 

Now, let us return to our method and calculate the principal operator of the same system in R, which is well defined 
and finite from the beginning of the problem. The result is 
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Hi 



dxdx'x^{x) / dt Kf{x,x')e 



dx dx' x^i^) 



dxi dx2 dx[ dx'2 



X (l)\x[) {x'^) / dt Kt{x[,x')Ktix' ,x'2) Ktixi,x) Kt{x,X2) e-*(^«-^V(xi) 0(x2) 



Xix') , 



+ 4/ dxidx2(^\xi) dtKt{xi,x')Kt{x',x)Kt{x,X2)e-*^"°-^^(j){x2) 



(127) 



where Kt{x,y) = ^71-7^171— • The condition (|102p gives 



(47rt)i/2 

J^^ dxdy %{x)My)KKx,y)e-'\''^^\ [ / da;' dy' ul{x')Kt{x' ,y')m,{y') 
(n-2)(n-3) 



+ 2(ri-2) / dt 



dt I dxdxidx2 ul{xi)ul{x2)Kt{xi,x) Kt{x2,x)il}o{x) 

2 

dxdyul{x)Kt{x,y)My) e-*l^-V / da;' dy' y')«o(y') 



n— 4 



Following the same analysis given above, we find the left hand side of (|128p for large values of \E, 



gr\ 



Ax |^o(x)|^ 



(8^t)l/2 \ 



dx' |«o(x')|' 



A 2^2^ 



(128) 



(129) 



'ar\ 



and the right hand side of it in the same limit, which is the analog of pi3p in one dimension, becomes less than the 
following term 



1 



Sl-E'grI (1 + 



nK[uo] 



dx \uo{x)\'^ 



2n 



\Egr\ (1^!!^^ 



(130) 



In one dimension, the Sobolev inequality for 2 < g < 00 is given as [45 



dx 



duo 



dx 



2 \ e 



dx luol 



l-Egrl 



> Si^ql / dx \uo\ 



2/q 



(131) 



where 0=i(l-f) 



and 



'(1-0)1- 



22/9(g- 2)(9-2)/'J 



q-2 ' 2 



(9-2)/? 



(132) 



with equality if and only if uo(a;) ~ ccosh ^^'"^ '^\b{x — a)) for some a G M, 6 > and c G C. Since we are looking 
for an upper bound to (|130p we will choose g = 4 so that = 1/4. Then the Sobolev inequality in (|13ip gives 



j dx \uo\^ < S^l(^J dx 



duo 



dx 



2 \ 1/2 



dx |wo|' 



3/2 



(133) 



where we have used the normalization of the wave functions and 6*1.4 = 3^/"*. Using this result in (|130p and from 
([T^ . we get 



< 



A-i/2[, 



2n 



A 2^2^- 2V3|£;3.|(1 + !^) |i?5.|(l + !^) 



(134) 



I-Egri 
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Keeping the leading order term on both sides, we obtain 

Let us define the variables z = nK[uQ\ and a = l/\Egr\, and then find the upper bound to the right hand side. This 
occurs at z = 1 /a so we get 

Egr > > (136) 

which is exactly the same result given in (jl26p in the leading order. Wc note that in this approach the kinetic energy 
of the center of mass motion is automatically set to be zero. We can also find the cigcnfunction from our analysis. As 
a result of the above theorem, the Sobolev inequality that we have used above is saturated if 



cosh(oj;) 

Here we have chosen the constant a = without loss of generality and the coefficient c = \fbj2 has been found from 
the normalization. The constant h can be determined from the solution z ~ nK\ai^ — \Egr\. Since the saturating 
solution ()137p satisfies 

dx|wo|^ = ^i^'/'M, (138) 

we obtain h — An/4. Therefore wc find exactly the same result obtained from the Hartree approximation (|125p . 
Incidentally, in this limit the wave functions could be taken as, 

uoix) = y^e-"^l"l/4 . (139) 
and they are related to the actual wave function of the system by our previous formula (j93p . 



/ 



V. RENORMALIZATION GROUP EQUATIONS 

The renormalization group equations (or Callan-Symanzik equations) for the system, where the particles do not 
interact with each other but interact with an external Dirac delta potential in two and three dimensional flat spaces, 
has been worked out in [l^, [13 ■ Many-body version of the same problem, where the particles interact via two-body 
delta potentials, has also been studied psl l26l. |49|. 

Recently, we have derived the generalization of the renormalization group equations of the above one-body model 
with N delta centers into two and three dimensional Riemannian manifolds |3l| . Here, we will show that the interacting 
version of the problem can be also studied explicitly, as we will see. 

One possible way for the renormalization scheme in order to determine how the coupling constant changes with the 
energy scale is to define the following renormalized coupling constant Xji{M) in terms of the bare coupling constant 
A(e) 



1 



Xr{M) A(e) 



dt- 



(140) 



where M is the renormalization scale (it is of dimension [i?]^/^). Then, the renormalized principal operator in terms 
of renormalized coupling constant is given by 



Hi 



Xr{M) 



d^xd^x' xlix) 

A12 JO 



dt 



Kf{x,x';g)e-*^"<^~^'> 



8TTt 



-5(2)(X,X') 



dlxdlx xlix) 



dlxidlx2dlx[dlx'2(l)l{x'-i)4il{x'2) / dt Kt{xi,x]g) Kt{x2,x\g) 



'0 

X Kt{x' ,x'^-g) Kt{x' ,x'2]g) e-'^"°-^^g{xi) (j)g{x2) + A / d^xi d^a;2 4>l{xi) 

oo 



X / dt Kt{x2,x;g)Kt{x',x;g)Kt{x',xi;g)e-'^"°-^Ug{x2) 



Xgix') . (141) 
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Here the bound state energies are again determined from the condition ^^{E)\^) = in the n-particle sector, 
however there is an ambiguity, we have a family of solutions for different choices of M and Xfi{M). To determine the 
value of Xii{M) at an arbitrary value of the renormalization point M, a natural choice would be to use the physically 

(2) 

measured two-body bound state energy Egr , if it exits, otherwise to use a scattering amplitude at some two particle 
energy. The solution then determines the relation between Ak(M) and M. Explicit dependence on M cancels the 
implicit dependence on M through Xji{M). In the case of two-body bound state energy, the principal operator acts 
on \0) (E) Jj^d^x ip{x)xl\0)- Hence, because of the condition for the bound states (|54p we obtain an equation, the 

(2) 

solution of which fixes Xji{M) as a function of M, Egr : 



1 



CO 



ij{x) = . (142) 



Even if we cannot explicitly solve this equation, the arbitrariness in the choice of the scale is reflected by expression 
below, 

^^d^HM,XuiM),E;g) ^ ^ 

dM ^ ' 

or 

^^7777 Afl(Af),ii;;5) = , (144) 



dM ' " ' BXr 



where 



(145) 

is called the /? function and equation (|144p is the renormalization group (RG) equation. This equation implies that 
the physics is independent of the choice of our renormalization scale. Using (|14ip in (|144p . we can find /3 function 
exactly 

/5(A«) = < • (146) 

This result is exactly the same as the one in fiat spaces given in the literature (26j so our problem is asymptotically 
free, too. 

We will now derive an analog of Callan-Symanzik equation for our principal operator and show that there is a 
simple solution of this equation, related to the fiow of the renormalized coupling constant. This will reconcile present 
method with the tools of conventional approach to field theories. 

In order to see this, we will use scaling property of the heat kernel in two dimensional Riemannian manifolds 

^t{x, y\ g) = j^^K^^2t{x, y; -f^^g) , (147) 

with the assumption that the manifold that we are interested in is stochastically complete, that is, d^x Kt{x, y; g) = 
1. There exists a unitary representation for the scaling transformation of the metric g i— ^~'^g such that the creation 
and annihilation operators transform like 



U(^)<j,g{x)UH^) = i-^^,-2g{x) , C/(7)0j(x)C/t(7) = 7"V^-2g(x) 

U{^)Xg{x)U\^) = ^-\^-.g{x) , U{^)xl{x)U\^)^^-\\_.^{x), (148) 

where we have used their commutation relations and the algebra of the orthofermions defined in (jlOp . Wave function 
normalization will be invariant under this transformation. 

Let us first simultaneously scale the energy by 7^ and the metric by 7~^ in the renormalized principal operator 
given explicitly in (|141|) and get 

^-{M,Xnm,l'E-,,--9)-^-^^^^^^^^^ 
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X / dt 

10 

1 



2 /_^/7- = S^d7-^9^ x\-2g{x) 



i^-2gXl d^-2gX2 d^-2gXl - 2 g ^2 _ 2 g ( 2: 1 ) 



X 0t-2„(a;2) / dtKtixux;-/ ^g) Kt{x2,x;-/ ^g)Ktix' ,x[;-f ^ g) Ktix' , x'^; j ^g) 



X e-*(^°-^'^V7-^9(^i) 07-^9(2^2) + 4 / d2_,^xi d2_,^a;2 0^_2g(a;i) 



X / dt Kt{x2,x- -,-^g) Ktix', x; -f-^g)Ktix', xi-n-^g)e-*^"''-''"''^ cj)^-2g{x2) 



X7 ^g{-^ ) 



Now we make a change of variable t y-^ j and use the scahng property of the heat kernel (|147p and obtain 



(149) 



dixdyxl-2gix) 



Kf{x,x';g)t 



-5f\x,x') 



2 /..O ^l^^l^' Xl-2gix) 



J^^ d^Xl dlx2 dlx[ dlx'^cl,\_._^{x[) 



^ <pl-2g(x'2) J dt Kt{xi,x;g)Kt{x2,x;g)Kt{x',x[;g)Kt{x',x'2;g) 



/•OO 

X / dii^t(a;2,a;;g)A't(a:',a:;5)ir*(x',xi;5)e-*^"'(^«-^'^V7-^s(a^2) 



dgXidgX2 (t)\-2g{Xi) 



X^-2g(a;') , (150) 



where we have used ^^,^-^2g(a^, a^') = I'^S^'' {x,x') and d; 



7 ^d^x. Using p48p . and inserting the identity 



U{"f)W{'^) in the appropriate places inside the above equation, we obtain for U'^{'^) $'^(M, A_r(M), 7^£'; 7 ^g) ?7 (7): 



t/t(7) ^^{M,\R{M),-i^E--i-^g) U{j) 



X I dt 
'0 



K^(i,i';g)e-"»'-':^ - 



Xg{x') 



d^Xid^a;2d^x'id^X20],(x'J(?!)t(x2) / dt Kt{xi, x; g) Kt{x2,x; g) 



X Kt{x\x[;g)Kt{x\x'^;g)e-'^"°-''^cl)g{xi)<l>g{x2) + A / d^xi d^xz 0],(xi) 



X / dtKt{x2,x;g) Ktix', x;g) Ktix', xi;g)e-'^""'^Ug{x2) 
Jo 



Xgix') 



where 



Therefore we finaUy obtain 



C/t(7)e 



-t7"'(ffo-7^ 



-t{Ho-E) 



C/t(7) cI>«(M, Afl(M),7'i?;7"'5) C/(7) = <i>^(7"'M, A«(Af),ii;;g) 



(151) 



(152) 



(153) 



It is important to note that we need to scale the metric as well. The idea of the metric scaling in deriving the 
rcnormalization group equation was motivated by [sof in the context of renormalization group in quantum field 
theory on curved spaces. Hence we have 



^d7 



C/t(7)$«(M, Afl(M),72i?;7-2,g)[/(7) = <J>«(7-iM, XniM),E;g) 
This leads to the renormalization group equation for W i'jj^-'^'iM, XRiM),j'^E;j~'^g)Uij) 



(154) 
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[/t(7)$«(Af, Afl(M),7'^;7"'5)C^(7) = 



If we postulate the following functional form for the principal matrix 

[/t(7)$«(A/,Ai^(A/),72£;;7-2g)t/(7) = /(7)$''(M, XRhM),E;g) 
and substitute into (|156p we obtain an ordinary differential equation for the function / 

d/(7) 



d7 







(156) 
(157) 
(158) 
(159) 



This has the solution 7(7) = 1 using the initial condition at 7 = 1. Therefore, we get 

t/t(7)$^(M, Afl(Af),72S;7-2.g)C/(7) = 'S^'^iM, Xr^M), E; g) , 

which means that there is no anomalous scaling. This interesting result has been derived in j25l [26j for the two-particle 
sector in flat space for T-matrix. 
By integrating 



I3{Xr) ^ M 



dXiiiM) _ XliM) 



dM 47r 

between M = M to M = we can find the flow equation for the coupling constant 

Xr{M) 



l + i^XR{M)\n^ 



(160) 



(161) 



Indeed, the above evolution can also be derived from the choice of our coupling constant given in ((T40| . One can 
explicitly check the relation (|159p if the coupling constant evolves according to (|16ip . First, we add and subtract a 
term in the time integral to {M, Xr(jM), E; g) (as indicated explicitly below) and use (|16ip : 



$^(M,Afl(7Af),£;;5) 

) 

dt 



Hi , Hi 



Xr{M) Att 
Kfix,x';g)e-'^""-^^ 



ln7- / d^a;d^a;'x^(a 



I M 



SlTt 



8TTt 



\i dgxdgx' Xg(a;) / d^ari d^xa dgx'^ d^^ 0g(a;i) '/'gl^) / dt Kt{xi,x]g) Kt{x2,x;g) 

^ J M'^ L J Jo 

X Kt{x',x[;g)Kt{x',x'2;g)e-"^"''-''Ua{xi)Mx2)+'i f d^xi d2a;2 (j^Uxi) 



we find 

'S>''iM,XRijM),E;g)-- 

} 

dt 



X / dtKt{x2,x;g)Kt{x',x;g)Kt{x',xi;g)e-*^"'>-''')^g{x2) 
d^xdlx'xlix) 



(162) 



Hi 



Ai?(Af) Jm 



Kf{x,x';g)e 



-t(Ha-E) 



8^ 



-42)(x,x') 



Xaix')^^ [ dlxdyxlix) 



dg.TidgX2dgxidgX2 ^t(x'i)(/)J(a;2) / dt Kt{xi,x]g) Kt{x2,x]g) 



X Kt{x\x'^-g)Kt{x',x'2]g)e-'^"''-''^g{xi)(t)g{x2)+'^ / d^xi dlx2 4>l{xi) 



X / dt Kt{x2,x;g) Kt{x' ,x;g) Kt{x' ,xi;g)e 



4>g{x2) 



Xg{x') 



(163) 



This is exactly equal to $^(7-11/, Xr{M),E; g) and this is indeed W{-f)<^^{M, Xr{M), j'^E; j^'^g)U{-f) due to (|T53| . 
This shows that one can alternatively find out evolution of the coupling constant which is given (|16ip from the scaling 
relation (|159p . 
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VI. CONCLUSION 



111 this paper, we have constructed a new non-perturbative renormahzation method to the many-body problem on 
two dimensional manifolds. The ground state energy is studied in the mean field approximation. The renormahzation 
group equation has been derived and the /3 function is exactly given, as a result it is shown that the model is 
asymptotically free. 
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